Tunable Kondo effect in double quantum dots coupled to ferromagnetic contacts 
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We investigate the effects induced by spin polarization in the contacts attached to a serial double quantum dot. 
The polarization generates effective magnetic fields and suppresses the Kondo effect in each dot. The super- 
exchange interaction (Jafm), tuned by the inter-dot tunnelling rate t, can be used to compensate the effective 
fields and restore the Kondo resonance when the contact polarizations are aligned. As a consequence, the direc- 
tion of the spin conductance can be controlled and even reversed using electrostatic gates alone. Furthermore, 
we study the associated two-impurity Kondo model and show that a ferromagnetic exchange coupling (Jpm) 
leads to an effective spin- 1 exchange-anisotropic Kondo model which exhibits a quantum phase transition in the 
presence of partially polarized contacts. 
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The study of spin-polarized transport in quantum dots 
(QDs), motivated by potential application for spintronic de- 
vices, has recently drawn a lot of attention both theoretically 
and experimentally J^l-Ql. Ferromagnetic contacts af- 
fect the dot charge dynamics so that spin-dependent tunnelling 
rates renormalize differently the dot level for each spin orien- 
tation HI 0. 0] . This is reflected in the appearance of an ef- 
fective magnetic field in the dot H-12] which suppresses the 
many -body Kondo state ll U llZl . routinely observed at low 



enough temperatures in QDs. However, this Kondo state can 
be restored by pro perl y tuning the QD level position (e) via 
a gate potential l^ llOll . Additionally, when the polarizations 
in the ferromagnetic electrodes are non-collinear, the induced 
effective magnetic field depends on the relative orientation of 
the two easy-axis, but it can, again, be compensated lfl3tl . 

Double quantum dots (DDs) have been extensively stud- 
ied because they offer interesting perspectives in quantum 



spintronics such as spin-based quantum computation 11411 . 
Pauli spin blockade 1 15 . 3, spin pumping lfl7ll . etc. More- 



over, DDs offer a natural experimental realization of the two- 
impurity Kondo problem (2IKP) 1 18-2(3]. The interplay be- 
tween the single dot Kondo effect and the interdot interaction 




24] and its ex- 
2611 . Coupling 
adds a further 



has been extensively studied theoretically 
perimental realization has been reported 
such DDs to ferromagnetic contacts (27 
experimental handle to the system, which is distinct from ap- 
plying an external magnetic field and may have interesting ap- 
plications in spin-dependent transport. 

In this letter we investigate the role of injecting spin- 
polarized charge carriers in serially coupled QDs (see upper 
inset in Fig. 1). We report pronounced transport phenom- 
ena due to the coupling with ferromagnetic electrodes, which 
are not present in single dots. In the absence of the particle- 
hole (p-h) symmetry, ferromagnetic contacts induce effective 
magnetic field in each QD. We show that for parallel magne- 



tization in the electrodes, these induced effective fields can be 
compensated by properly tuning the inter-dot tunneling am- 
plitude t, thereby restoring the Kondo states in each QD. Fur- 
thermore, we demonstrate the behaviour of the DD set-up as 
a controllable spin-filter device by gating the dots and varying 
either the lead polarization or, which is more useful, the inter- 
dot tunnelling t. Finally, we complete our study by analysing 
the role of the polarization on the 2IKP. In particular, we find 
that the equivalent model with a ferromagnetic spin-spin inter- 
action Jfm behaves as a spin-1 exchange-anisotropic Kondo 
model and exhibits a Kosterlitz-Thouless quantum phase tran- 
sition at finite polarization. 

Model - The double dot set-up (Fig. 1) is modelled as a 
two-site Hubbard model. The Hamiltonian reads: 

U = E £akacl ka C ak a + E £«a<r + U n a ^n a ^ 



k.a 



+ t 



E { d l d ^ + h - c -) + E (^4<a- + h - c -) > (d 



Here c Q fe CT annihilates an electron in the electrode a e {1,2} 
with wave-vector k and spin a =f , L Similarly, d aa destroys 
an electron with spin a on the dot a. Each dot is also con- 
nected to a contact with hybridization amplitude V a . Ferro- 
magnetic contacts are described by employing spin-dependent 



tunnelling rates: T Q . 



n\V a \ v ai j, v aa being the spin- 



dependent density of states at in the contacts. The polar- 
izations are parametrized by p a , defined through T a ^ = 
T(l + p a )/2, T„4 = T(l - p a )/2; we assume the total hy- 
bridization strength T to be constant and equal for both con- 
tacts. In this work we will consider the case of collinear polar- 
ization between the leads with parallel (P) (pi = P2 = p) and 
anti-parallel (AP) alignment (pi = — p2 = p). We solve H 
using the numerical renormalization group technique (NRG) 
& 

Particle-hole symmetric case, 5 = e + U /2 = 0. - The 
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Figure 1. (Color online) Upper inset: Illustration of a serially 
coupled double dot set-up for the parallel P configuration [major- 
ity(minority) spins of both reservoirs are down(up) spins, t is 
the inter-dot tunneling amplitude, V is the lead-dot hybridization 
strength. From (a) to (f): up and down spectral densities for the dot 1. 
Parameters are e = —U/2, U = 7, and T — 1. (g) Spin-spin corre- 
lation function for the P arrangement (1 and 2 sub-indices denote the 
left and right dot, respectively), (h) Double dot phase diagram for the 
P case when t and p are varied: the three phases are the Kondo phase, 
the anti-ferromagnetic phase (AFM) and the molecular-orbital phase 
(MO) in which for a very large t the two dots behave as an effec- 
tive one-dot system with energy levels at e ± t. In this case, the usual 
spin-1/2 Kondo physics emerges. Lower inset: Double dot system in 
the AP configuration [majority (minority) spins in the left reservoirs 
are down(up) spins and majority(minority) spins in the right reser- 
voir are spin-up(down) electrons]. Up and down spectral densities 
for the dot 1 in the AP lead orientation are plotted in (i), (j), (k), and 
(1) . Parameters are e = -U/2, (7 = 7, and T = 1. 



dot-1 spectral functions Ai a (u) for the P and AP alignments 
are shown in Fig. 1. For p = and small t [Fig. 1(a)], 
the spectral functions for both spin alignments show a single 
peak (the Kondo resonance) pinned at the Fermi level. As t in- 
creases, one moves from the regime dominated by the Kondo 
effect to a phase governed by the antiferromagnetic (AFM) 
coupling between the spins [Figs. 1(d), (e), and (f)]. The 
AFM regime is evidenced by a double-peak DOS with peaks 
at oj ps ±J/2 where J = Jafm = 4i 2 /(7 is the superex- 



change interaction. For finite polarization, p ^ 0, each dot is 
predominantly affected by the polarization of the neighboring 
contact and only indirectly (thus weakly) by the other contact, 
therefore P and AP arrangements show similar behavior [cf. 
Fig. l(a)-(f) and Fig. l(i)-(l)]. For small t, in the P orien- 
tation, [Fig. 1(c)] the spectral weight at Ep becomes spin- 
dependent, i.e., Aif(Ep) ^ Aii(Ep) but there is no spin- 
splitting. The width of the Kondo peak reduces according to 
[ll T K (p) a DexpU-l/foJK + v± J K )][twh- X {p) / p]} 
where vJk = ST/ ttII at the p-h symmetry point. Thus, by 
increasing p [Fig. 1(e)], the Kondo temperature is lowered. 
When Jafm ^ 27^ (p), the Kondo effect is destroyed and the 
dot spins bind into a local singlet which is reflected in a split 
spectral density [Fig. 1(e)]. The correlation function (Si • S2) 
which measures the dot spins alignment is shown in Fig. 1(g) 
for increasing p for t = —0.769. Accordingly, the spectral 
function exhibits a splitting of magnitude Jafm [Fig. 1(e)] 
signaling the formation of the AFM singlet state. This is the 
same type of cross-over that one observes in the p = model 
as a function of t. This result implies that there is a continuous 
cross-over line in the (p, t) plane separating the Kondo phase 
from the AFM one. A schematic phase diagram summarizing 
this result is depicted in Fig. 1(h). Spectral dot densities in 
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Figure 2. (Color online) Spectral function of dot 1 when the system is 
away from the particle-hole symmetric point. Parameters are [7 = 7, 

r = 1, e = -U/2 + 1. 

the AP orientation are shown in Figs. l(i)-(l). As mentioned 
previously, they exhibit a similar behaviour than in the P con- 
figuration. We find remarkable differences between the two 
magnetic orientations (P) and (AP) in the half-metallic limit- 
ing case when p — > 1 [illustrated in Fig. 1(e), and (f) for the 
P case and Fig. l(k), and (1) for the AP case]. For the P case 
the spin-down spectral density at the Fermi level Ep (ut = 
in the plots) vanishes, while for AP it remains finite. Con- 
versely, for the P arrangement the spin-up spectral function is 
finite at Ep, while for AP it goes to zero. The difference is 
due to an interference effect that can be studied analytically 
in the non-interacting ((7 = 0) limit, with the result holding 
more generally. For parallel alignment (p — > 1), spin-down 



3 



electrons are localized inside the DD system, while the spin- 
up electrons can still flow between the contacts; the U = 
spectral function for spin-down electrons has delta-like peaks 
at cj = ±t and is zero elsewhere, while the spin-up spectral 
function is finite at Ep. For AP, in the p = 1 limit the spin-up 
electrons from the contact 1 can enter both dots, but cannot 
exit at the right contact. The U = spin-up spectral functions 
thus shows a peak at Ep for the dot 2 and, consequently, a null 
spectral density at Ep for the dot 1 due to interference. For 
the other electron spin, the behavior in the two dots is simply 
reversed. The difference in the p —> 1 limit becomes even 
more pronounced as t increases. 

General case, j ^ 0. - Away from the p-h symmet- 
ric point, the spin splitting is generated by virtual processes 
proportional to the spin-dependent hybridization functions: 
B c s = Sef — 6e± (with a = —a) 

^^-l/^/^-^l+^p)). (2) 

7T J \ id — 6 e+U — LO ) 

Accordingly, at p ^ and for general values of parameters, 
the spectral function exhibits spin splitting (see Fig. 2). For 
small t, this splitting is fully analogous to that observed in 
single QDs [e.g., see Fig. 2(a), and (b)]. For a single QD, 
in the AP case there is no induced field because of the di- 
rect compensation of the contributions from both contacts. In 
the DD case, however, there is an induced field in each dot 
for both P and AP cases, but they differ in the direction of 
the fields: they are aligned along the same direction for the 
P case (-Bi iC ff = i?2,cft), while they point in opposite direc- 
tions for the AP case (Bi c s = — i?2,eff)- In single QDs, the 
splitting can only be restored by the application of an external 
magnetic field. In DDs, we find, in contrast, that the splitting 
compensation can also be achieved by the exchange fields due 
to the inter-dot exchange coupling Jafm: for a specific value 
of t the compensation occurs [Fig. 2(c)] and beyond this t 
value the splitting is shown again [Fig. 2(d)]. This can be 
understood as follows: seen from dot 1(2), the exchange cou- 
pling Jafm between the dots can be regarded as an effec- 
tive magnetic field: -Boxc = </afm<S Z 2(i). The restoration is, 
however, only possible in the P case while in the AP case the 
splitting only grows when t is increased [see Fig. 2(e) and 
(f)]. Note that this type of restoration of the Kondo peak has 
also been predicted in the two-impurity Kondo problem in an 
external magnetic field |3^, 31] and verified experimentally 
1 3211 . Interestingly, all these features shown in the dot spec- 
tral density are reflected in a measurable transport magnitude: 
the linear conductance. We compute the conductance through 
the DD system as a function of the inter-dot tunnelling t and 
the dot level position e for three values of p in Fig. 3 (upper 
panel). For p = we observe the standard results for the 
conductance of a DD: along the particle-hole symmetric line 
the conductance is low in the small-i and large-i limits, and 
it peaks at the cross-over point Jafm ~ 2Tjf. For very large 
t, the conductance becomes high when one of the molecular 
orbitals at energies e ± t is tuned to the Fermi level signaling 
the molecular orbital phase. The cross-over conductance peak 
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Figure 3. (Color online) Conductance plots versus 5 = e + U/2 and 
t. Upper panel: Linear conductance for p = 0, 0.6, 0.9 in the par- 
allel arrangement. Lower panel: spin-up (left), spin-down (center) 
conductance and spin conductance (right) for p = 0.6. Parameters 
are U = 7, T = 1. 



and the molecular-orbital conductance peak are smoothly con- 
nected in the (t, e) plane by a high-conductance ridge which 
reaches the unitary conductance, G = 2e 2 /h. The influence 
of polarized contacts is more dramatic for the P alignment, 
due to the interplay between Jafm and the effective fields. 
The DD can namely behave as a spin-filter device for a a fi- 
nite polarization in the P arrangement. This can be seen from 
the spin-up conductance Gf, the spin-down conductance Gj,, 
and the spin conductance G s = G-\ — Gj,, which are shown in 
the lower panel of Fig. 3. For finite p the conductance ridge 
no longer reaches the unitary limit. The maximum conduc- 
tance coincides with the restoration of the Kondo effect when 
B cxc + B c g = 0. Gi exhibits a sharp peak in the (t, e) plane; 
the narrow width is due to the weak hybridisation of the spin- 
down electrons that reduces the Kondo scale as p increases. 
As a consequence, one finds that the direction and amplitude 
of the spin conductance G s can be tuned by the parameters e 
and t (shown in Fig. 3). Therefore the DD acts as an efficient 
spin-filter device with potential applications in the construc- 
tion of spintronic devices. 

Equivalent two-impurity Kondo model. - We also study 
the associated Kondo model with direct inter-dot magnetic ex- 
change coupling of the form % = JSi • S2 and no hopping 
term, i.e., the hopping term t is replaced by a pure exchange 
term with J = Jafm = 4i 2 /£/. Provided that the charge 
transfer term between the leads (which corresponds to a co- 
tunneling term through the two dots) can be neglected, one 
finds a true quantum phase transition which corresponds to 
the critical behavior in the 2IKM 12411 . We found that polar- 
ization in the leads does not affect this conclusion. However, 
in an experimentally realistic system, this phase transition is 
replaced by a smooth cross-over as a function of t 13311 . 

For completeness, we have also analysed the case of a ferro- 
magnetic exchange interaction, i.e., J = — Jfm < 0. The cor- 
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responding spectral functions are shown in Fig. 4. For normal 
electrodes and small Jfm, each dot exhibits a single-channel 
spin-1/2 Kondo state, while for a large Jfm> the dot spins 
first rigidly bind at some high temperature into a spin-1 state 
which is then screened at a lower temperature. Note that the 
shape of the Kondo resonance is therefore different for small 
and large- Jfm cases [corresponding to Fig. 4(a) and Fig. 4(b), 
respectively]. At finite p [Fig. 4(c), and Fig. 4(f)], the Kondo 
effect is destroyed and the NRG results show that the DD set- 
up has a residual In 2 entropy, indicating the occurrence of a 
quantum phase transition from the Kondo screened phase to 
a partially-screened phase. The transition is of the Kosterlitz- 
Thouless type and it occurs at a value of p which depends on 
Jfm- Thus, there is a line of phase transitions in the (p, Jfm) 
parameter plane. This behavior is not unexpected. It should 
be recalled that for a single QD, the presence of spin polariza- 
tion results in the exchange-anisotropy of the effective low- 
energy Kondo model, J± ^ J z . The same effect is expected 
in the present spin-1 Kondo model. It is known, however, that 
anisotropic high-spin Kondo model feature partially-screened 
states in their phase diagrams 0411 . 
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Figure 4. (Color online) Spectral functions for the model with ferro- 
magnetic exchange coupling for three different polarizations p — 0, 
p = 0.6 and p — 0.9 (P case) and two distinct Jfm strength cou- 
plings. 

Conclusion. - Ferromagnetic contacts alter profoundly the 
transport properties of serially coupled double dots. The inter- 
play between the anti-ferromagnetic exchange coupling and 
the effective fields induced by the polarized leads in the paral- 
lel contact configuration may reinforce the Kondo effect. By 
tuning the interdot tunnelling and the dot gates the spin con- 
ductance reverses its sign, thus the DD acts as a spin-filter 
device. We have also shown that the 2IKM with ferromag- 
netic interaction exhibits a Kosterlitz-Thouless type transition 
ascribed to the partially spin polarized contacts. We propose 
double dot carbon nanotubes attached to polarized contacts as 
the best candidate to observe our predictions. In carbon nan- 
otube quantum dots, large polarizations and much stronger 
Kondo states have been observed in comparison with semi- 
conductor quantum dots IB5I438I1 
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I. EQUATION OF MOTION STUDY OF THE NON-INTERACTING CASE, U = 

We consider the [7 = limit of the Hamiltonian presented in the main text. Our goal is to compute the spectral functions 
using the equation-of-motion method: 

z((A-B)) z = ([A,B]+) + (([A,H]-;B)) Z (1) 

where z is the frequency parameter in the complex plane, ((A; B)) z denotes the correlator between the operators A and B, 
while [A, B]- and [A, B} + are a commutator and an anticommutator, respectively. We introduce the notation Gi a .j T (z) = 
((dio-; dj T )) z for the impurity Green's functions, where i,j <E {l,2}andcr,T 6 {ti^}- Since the spin components along the 
z-axis are conserved, only a = r parts of the Green's function are non-zero, thus we may also write Gij i<T (z) = Gi tT j (J (z). 

Furthermore, we define G a ki,a(z) = ((c a ka', d\ a )) z where k is a wave-vector in the lead a S {1, 2}. 
We obtain (omitting the argument z of Green's functions for brevity) 



(2) 



and 



From the latter equation it follows 



zGn,c, — 1 + eGii.c + tGii^ a + V\ Giki,o, 

fc 

zG\2 t a = tG\2 t a + tG22,<r + V\ G\k2,cn 

k 

zG21,a = tG2l,a + tGn^ + V% G^fcl.crj 

k 

zG22,cr = 1 + eG22,cr + tGl2,cr + Vi G2k2,o, 



zG otki,(T — £-otk<jG ctki,cr ~t~ VqG ai.a ■ (3) 



Gaki.a — G a i j(J . (4) 

Z taka 



We introduce the (spin-dependent) hybridization function as 

V 2 

A ^ = E _° ■ ^ 

k 

The spin dependence arises from the er-dependence of the band dispersion e a ka-, i.e., from the spin dependence of the conduction- 
band density of states. 

Solving the system of equation we obtain 



G22,a{z) 



z — e — Aicr -~r — 

1 

z — e — A 2 .r — 



(6) 



The spectral functions are then simply obtained as 

A i<7 (u) = — ImG ii)ff (w + id), (7) 

where S — > + . 

We now consider the half-metallic limit, p4l. For parallel alignment, this implies 

Ai t = A 2t = A > A H = A n = °- (8) 
We take the wide-band limit where A = — iT. We obtain 

p ^2 _|_ j2 _|_ p2 

A lt («) = MM - - [(w _ e + t)a + ra][(w _ e _ t)a + ra] . W 
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and 



A x4 .(w) = = i<J(w - e - t) + ^5(cj -e + t). (10) 

For antiparallel alignment we have 

A lt = A 2i = A, A u = A 2t = 0. (11) 



It follows 



and 



A n (u}) = A 2i (uj) = - ^ ^ -^r (12) 

7r i 4 + 2t i {uj — e) z + (u — e) [i + (to — t) \ 



Note that at the particle-hole symmetric point (e = 0), the spectral functions Ai^ has a zero while A 2 ^ is finite at u> = 0. This 
result also holds in the interacting case with U ^ 0. 



